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CHAPTER  I 
INTRODUCTION 


This  report  grew  out  of  work  at  Applied  Research  Laboratories, 
The  University  of  Texas  at  Austin  (ARL:UT),  In  ocean  acoustics, 
specifically  acoustic  propagation  between  source  and  receiver  located  in 
the  water.  In  shallow  water  regions  this  propagation  Involves  a  sig¬ 
nificant  amount  of  Interaction  with  the  seafloor.  This  Interaction 
affects  the  amount  of  loss  of  signal  level  incurred  as  a  signal 
propagates  from  source  to  receiver.  This  report  considers  the  acoustic 
attenuation  In  the  seafloor  sediments  and  the  effects  this  attenuation 
has  on  acoustic  propagation  In  shallow  water.  Also  presented  Is  a 
preliminary  technique  for  estimating  these  attenuations  as  a  function  of 
both  acoustic  frequency  and  sediment  depth.  This  technique  does  not 
require  physical  sampling  of  the  seafloor,  but  uses  acoustic  field 
measurements. 

The  estimates  of  the  attenuation  In  the  sediments  become  part 
of  a  geoacoustic  description^  of  the  measurement  region.  These 
geoacoustic  descriptions,  comprised  of  sound  speed,  attenuation,  and 
density,  together  with  water  mass  parameters,  are  used  to  numerically 
compute  acoustic  propagation  In  the  water.  These  calculations  can 
provide  propagation  loss  for  source  depths,  receiver  depths,  and 
freipiencles  not  covered  In  the  measurements,  thus  providing  a  complete 
acoustical  analysis  of  the  measurement  region. 

A.  Shallow  Water 


In  ocean  acoustics  the  term  "shallow  water"  has  been  given 
different  definitions  by  different  authors.  Brekhovsklkh^  and  Tolstoy 
and  Clay*  define  shallow  water  In  terms  of  the  wavelength  of  sound  at  the 
frequencies  of  Interest.  Other  authors,  such  as  Officer,^  define 


shallow  water  by  a  physical  water  depth.  The  definition  of  shallow  water 
used  In  this  report  was  obtained  by  first  considering  the  frequencies  of 
Interest  In  related  research  at  ARL:UT^*^  and  then  choosing  a  specific 
depth  as  the  maximum  depth  to  be  considered.  The  frequencies  of  Interest 
range  from  a  few  tens  of  hertz  up  to  2000  Hz.  Within  this  frequency 
range  the  compresslonal  wave  attenuation  Is  the  dominant  loss  mechanism 
In  the  sediment.  For  this  application,  the  definition  of  shallow  water 
has  been  chosen  based  on  the  wavelengths  of  the  lower  frequencies.  The 
characteristic  depth  was  chosen  to  be  500  m  and  shallow  water  Is 
Identified  as  those  regions  of  water  depth  less  than  500  m. 

This  shallow  water  definition  Includes  not  only  the  continental 
shelves,  but  also  the  shallow  seas  which  border  the  continental  margins. 
These  regions  are  regions  of  variability.  Major  variations  occur  In  the 
water  depth,  water  mass,  and  bottom  composition.  These  variations  will 
occur  within  a  region,  as  well  as  between  regions.  The  variations  In  the 
environment  result  In  variations  In  the  acoustic  propagation.^  Seasonal 
changes  In  the  sound  speed  profile,  changes  or  differences  In  water 
depth,  and  differences  In  bottom  composition  all  affect  the  propagation 
by  modifying  the  bottom  Interactions.  When  the  seasonal  changes  result 
in  Increased  bottom  Interaction,  the  propagation  losses  Increase. 

In  shallow  water  where  the  bottom  has  a  significant  Impact  on 
the  propagation,  variations  In  the  sediment  parameters  should  not  be 
Ignored  In  an  analysis  of  the  region.  In  particular,  the  variations  In 
attenuation  within  the  first  50  m  of  sediment  can  affect  propagation. 
Previous  authors  have  Ignored  depth  variability  In  attenuation.  Those 
authors  concluded  that  the  attenuation  was  not  linearly  related  to 
frequency.  They  also  concluded  that  shear  wave  excitation  was  an 
Important  loss  mechanism;  however,  shear  wave  speeds  significantly 
greater  than  those  measured  In  the  sediments  were  required  In  their 
analyses.  Ignoring  the  depth  variability  of  attenuation  can  therefore 
result  In  Incorrect  Interpretation  of  the  frequency  dependence  of  the 
attenuation. 


The  Idea  of  developing  a  procedure  for  estimating  the 
attenuation  profile  In  shallow  water  regions  grew  out  of  the  analysis  of 
acoustic  propagation  In  deep  water. Methods  for  measuring  the 
attenuation  profile  In  deep  ocean  sediments  were  developed  to  aid 
analysis  of  acoustic  propagation  data.  The  resulting  attenuation 
profiles  were  then  used  In  the  acoustic  description  of  the  sediments  In 
the  measurement  regions.  As  the  remainder  of  this  report  will  show.  It 
Is  necessary  to  use  an  attenuation  profile  In  the  acoustic  description  of 
the  sediment.  Therefore  a  method  was  needed  to  estimate  the  attenuation 
profiles  In  these  shallow  water  regions  using  the  measured  propagation 
data. 

g 

Hamilton  has  compiled  results  from  acoustic  measurements  which 
show  that  attenuation  does  vary  with  depth  Into  the  sediment.  A 
significant  contribution  to  this  compilation  was  based  on  work  presented 
by  Neprochnov,^^  who  was  able  to  analyze  deep  ocean  propagation  data  In 
terms  of  the  depth  variations  In  the  acoustic  properties.  Following 
Hamilton,  Mitchell  and  Focke^  and  later  Jacobson  et  a1.^^  have  presented 
measured  attenuation  profiles  In  deep  ocean  sediments. 

Neprochnov,  unfortunately,  did  not  present  the  method  used  to 
obtain  his  results.  However,  both  Mitchell  and  Focke  and  Jacobson  et  a1. 
did  discuss  the  methods  Involved.  Although  these  two  procedures 
differed,  both  techniques  were  based  on  multipath  analysis.  The  received 
signal  from  an  Impulsive  source  can  be  time  separated  Into  multipath 
arrivals,  each  with  Identifiable  ray  path  histories.  Figure  I.l  Is  an 

O 

example  of  a  received  signal  recorded  In  deep  water.  Also  presented  In 
this  figure  are  the  predicted  ray  arrivals  and  the  associated  ray  path 
histories.  The  histories  Include  the  number  of  bottom  bounces,  the 
bottom  grazing  angle,  and  a  reference  loss.  Measured  values  presented 
In  the  figure  Include  the  measured  loss  and  the  loss  per  bottom  bounce 
(bottom  loss)  for  each  arrival. 


The  bottom  loss  values  contain  information  about  the 
attenuation  profile  in  the  sediment.  Given  a  sound  speed  profile  in  the 
sediment,  ray  paths  through  the  sediment  can  be  identified.  The  bottom 
losses  then  provide  estimates  of  the  energy  lost  along  each  path,  and 
path  lengths  obtained  from  ray  tracing  yield  estimates  of  the  attenuation 
profile. 


These  deep  water  methods  can  be  seen  to  depend  upon  the  time 
separability  of  the  received  signals.  To  apply  these  procedures  to 
shallow  water  data  would  require  a  time  separability  in  the  data. 
Figure  1.2  presents  a  signal  received  in  a  shallow  water  region  with  a 
clay  sediment.  At  50  Hz  the  signal  is  composed  of  overlapping  arrivals 
and  the  total  received  signal  spans  a  time  of  approximately  4.5  sec.  At 
200  Hz,  the  signal  presents  a  different  structure;  there  is  one 
significant  peak  and  the  total  duration  is  approximately  2.0  sec. 

Figure  1.3  presents  a  received  signal  in  shallow  water  which 
has  propagated  over  a  sand  bottom.  The  source  waveform  is  composed  of  an 
initial  pressure  impulse  followed  by  three  secondary  impulses  at 
intervals  of  approximately  40  msec.  At  both  100  and  400  Hz,  the  received 
signal  is  made  up  of  a  single  arrival  with  this  pulsed  waveform. 

In  both  shallow  water  cases  above,  the  multipath  structure 
required  by  the  deep  ocean  attenuation  analysis  procedures  is  not 
present.  General  ray  theory  representation  of  propagation  used  in  those 
procedures  cannot  describe  the  characteristics  seen  in  the  data.  To 
estimate  the  attenuation  from  propagation  data,  the  propagation  needs  to 
be  analyzed  by  a  wave  theory  model,  such  as  a  normal  mode  model.  The 
problem  in  estimating  the  profile  will  be  associating  measured 
propagation  characteristics  with  depth  dependent  mode  characteristics. 

C.  Sunwary  of  the  Report 


The  remainder  of  this  report  is  divided  into  three  chapters. 
Chapter  II  is  on  acoustic  attenuation  in  the  sediments.  Chapter  III 
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FIGURE 

ENVELOPE  OF  SHOT  IN  V3  OCTAVE  BANDS 
PROPAGATED  OVER  A  CLAY  BOTTOM 
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FIGURE  1.3 

ENVELOPE  OF  SHOT  IN  V3  OCTAVE  BANDS 
PROPAGATED  OVER  A  SAND  BOTTOM 
RANGE  -  158  km,  SOURCE  DEPTH  -  91  m, 
RECEIVER  DEPTH -110m 


ARL:UT 


presents  the  preliminary  technique  for  estimating  the  attenuation 
profile.  Finally,  attenuation  profiles  are  estimated  for  three 
measurement  regions  In  Chapter  IV. 

The  discussion  on  attenuation  In  Chapter  II  notes  that  there 
are  two  points  of  view  on  the  variability  of  the  attenuation.  The  point 
of  view  taken  In  this  report  Is  based  on  measured  results:  the 
attenuation  varies  linearly  with  frequency  and  It  varies  with  depth.  The 
method  for  Including  the  sediment  attenuation  In  the  calculation  of 
propagation  loss  Is  also  presented.  Calculations  are  then  used  to 
demonstrate  how  a  nonlinear  frequency  dependence  of  attenuation  can  be 
obtained  If  depth  variations  In  the  attenuation  are  Ignored. 

The  technique  for  estimating  the  attenuation  profile  as 
presented  In  Chapter  III  Is  based  on  the  normal  mode  theory  and  uses 
propagation  loss  measured  at  several  frequencies.  Shallow  water 
propagation  In  the  frequency  range  below  approximately  2  kHz,  beyond  a 
few  water  depths.  Involves  only  the  low  order  normal  modes.  For  these 
modes,  there  Is  a  correlation  between  the  acoustic  frequency  of  the 
propagating  signal  and  the  amount  of  penetration  Into  the  sediment.  At 
the  high  frequencies,  the  modes  penetrate  only  the  surface  layer.  As  the 
frequency  decreases,  the  modes  penetrate  to  the  deeper  layers.  By 
calculating  the  normal  modes  and  determining  a  depth  of  penetration  as  a 
function  of  frequency,  attenuation  values  are  estimated  as  a  function  of 
depth. 

In  Chapter  IV  attenuation  profiles  are  estimated  for  three 

shallow  water  regions.  In  each  region  the  attenuation  values  were  within 

the  range  specified  by  Hamilton  for  the  sediment  types  found  In  the 

measurement  regions.  In  the  sand  sediment,  the  attenuation  was  found  to 

decrease  with  depth.  Again,  this  result  was  In  agreement  with  laboratory 
12  13 

results  and  theoretical  results.  In  the  clay  regions,  one  profile 
was  found  to  Increase  with  depth  while  the  other  profile  decreased  with 

Q 

depth.  Based  on  a  compilation  of  measured  results,  the  attenuation  In 
homogeneous  clay  sediments  generally  Is  expected  to  increase  with  depth. 


The  differences  between  the  two  estimated  profiles  and  the  generally 
expected  results  were,  however.  In  agreement  with  large  variability  of 
attenuation  In  these  sediment  types. 

These  linear  frequency  dependent  compress lonal  attenuation 
profiles  were  sufficient  to  account  for  the  frequency  characteristics 
seen  In  the  measured  propagation  loss  results.  Additional  loss 
mechanisms,  such  as  shear  waves  In  the  sediments,  were  not  required. 
Accurate  modeling  of  the  acoustic  propagation  and  an  understanding  of  the 
acoustic  Interactions  with  the  bottom  resulted  in  attenuation  profiles 
which  agreed  with  previously  reported  results.^ 


CHAPTER  II 

OCEAN  SEDIMENT  ATTENUATIONS 


The  conpresslonal  wave  attenuation  due  to  the  seafloor 
sediMnts  can  have  a  significant  Influence  on  acoustic  propagation  In 
shallow  water.  At  low  frequencies,  where  the  sediment  Is  the  only 
significant  source  of  attenuation,  this  parameter  determines  how  rapidly 
propagating  acoustic  signals  decRy  with  range.  Before  discussing  a 
method  of  estimating  these  attenuations,  the  attenuation  In  the  sediment 
should  be  discussed.  This  will  provide  an  understanding  of  the 
characteristics  of  the  attenuation  Including  the  range  of  values  to  be 
eapected.  There  will  also  be  a  discussion  on  how  attenuation  enters  Into 
the  propagation  loss  calculations  since  these  calculations  are  often  used 
to  ehtain  estimates  of  attenuation.  Thore  will  also  be  a  discussion  on 
the  effects  of  Ignoring  certain  characteristics  of  the  attenuation  In 
terms  of  data  analysis  results. 

A.  Proeartlos  of  Sediment  Attenuation 

Acoustic  compresslonal  wove  attenuation  within  a  sediment 
varies  with  acoustic  frequency,  with  depth  Into  sediment,  and  with 
sediment  type.^  How  the  characteristics  of  attenuation  are  related  to 
each  of  these  parameters  Is  of  Interest  In  estimating  the  attenuation 
profile  since  the  Interplay  between  these  relationships  contribute  to  the 
frequency  depe  ice  of  the  pr<H>agat1on.  The  effects  on  propagation 
resulting  frc  attenuation  varying  with  depth  and  frequency  are 
Interrelated,  estimate  of  one  relationship  Is  dependent  upon 
assumptions  ab  ^  the  other.  A  knowledge  of  the  variations  with  sediment 
type  are  of  Interest  when  estimates  are  obtained;  that  Is,  If  the 
estimates  are  to  be  accepted,  they  need  to  be  consistent  with  values 
eapected  for  the  sediment  type.  How  the  attenuation  actually  varies  with 


each  of  these  parameters  has  not  been  fully  established.  There  are 
basically  two  points  of  view:  one  based  on  In  situ  measurements,  and  the 
other  based  on  theoretical  considerations.  They  differ  In  some  very 
Important  ways. 

The  relationship  between  the  acoustic  frequency  and  the 
sediment  attenuation  Is  one  point  of  controversy.  The  mathematical 
models  based  on  the  theoretical  work  of  Biot,^^  Stoll  and  Bryan, and 
Stoll, consider  two  loss  mechanisms:  frictional  losses  between  the 
grains  and  fluid  viscous  losses  of  the  pore  fluid.  The  model  predicts  a 
linear  relation  between  frequency  and  frictional  attenuation.  For  the 
viscous  losses,  however,  the  model  predicts  a  relationship  that  changes 
from  f  at  low  frequency  to  f  '  at  high  frequency.  In  the  frequency 
range  of  Interest,  the  actual  relationship  between  the  total  attenuation 
and  frequency  will  be  a  combination  of  these  two  mechanisms  and  Is 
dependent  upon  the  sediment  type. 

For  sediments  of  high  porosity  (relative  volume  of  the  fluid, 

given  In  percentages),  such  as  clays  and  silts,  the  frequency-attenuation 

relationship  Is  almost  linear.  For  low  porosity  sediments,  such  as 

sands,  the  viscous  losses  dominate  the  attenuation  and  the  model  predicts 
2 

an  f  relationship  In  the  low  end  of  our  frequency  range,  changing  to  an 
f^/2  relationship  at  high  frequencies.  The  specific  relationship  at  any 
selected  frequency  Is  dependent  upon  the  grain  size. 

Interpretation  of  experimental  acoustic  propagation  measure¬ 
ments  have  Indicated  nonlinear  relationships  between  frequency  and 
attenuation  given  as 


a»Kf' 


Based  on  propagation  loss  measurements,  Rubano^”  found  the  attenuation  to 
be  proportional  to  the  1.1  power  of  frequency  (n»l.l).  Ingenlto,^^ 


also  working  with  propagation  loss  data,  found  attenuation  related  to 
the  1.75  power  of  frequency  (nal.75).  Analyses  of  acoustic  propagation 
data,  using  the  Blot-Stoll  model  to  provide  the  sediment  attenuation, 

1  o 

have  resulted  In  similar  relationships.  The  analyses  leading  up  to 
the  nonlinear  relationships  have  not  considered  depth  variations  In  the 
attenuation.  As  will  be  shown  later,  this  oversight  can  lead  to 
erroneous  Interpretations. 
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The  theoretical  work  of  Blot  assumed  a  uniform  spherical 

IQ 

pore  size.  Hovem*^^  has  shown  that  the  attenuation-frequency 
relationship,  due  to  viscous  losses.  Is  dependent  upon  this  assumption. 
Using  a  modified  suspension  model,  where  calculations  are  made  as  the 
grains  were  In  suspension  and  there  Is  no  effect  of  Intergrain  contact, 
Hovem  showed  that  for  mixed  grain  sizes  the  attenuation  tends  toward  a 
linear  relationship  with  frequency  over  a  broad  frequency  range,  while 
still  agreeing  with  the  Blot  model  for  uniform  grain  sizes. 

Hamilton  has  compiled  and  analyzed  as  many  as  possible  of  the 
attenuation  measurements  taken  In  unllthifled  marine  sediments.^  He 
found  that  the  attenuation  Is  approximately  linear  with  frequency. 
Unfortunately  there  Is  a  paucity  of  data  In  the  frequency  range  10- 
1000  Hz,  and  one  may  need  to  rely  upon  the  modeling  to  fill  the  gaps. 
Fortunately,  the  result  of  Hovem  Indicating  that  a  linear  relationship 
can  be  assumed  across  our  entire  frequency  range  Is  In  agreement  with 
the  data. 

At  any  given  frequency,  the  attenuation  In  different  sediments 
Is  found  to  vary  by  at  least  an  order  of  magnitude.  Assuming  a  linear 
frequency  dependence,  Hamilton  was  able  to  exhibit  an  attenuation 
dependence  related  to  sediment  type.  The  general  trend  Is  for  the 
attenuation  to  Increase  with  decreasing  porosity.  It  ranges  from  0.01- 
0.10  dB/kHz-m  for  clays  up  to  0.2-0. 9  dB/kHz-m  for  sands. 
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In  addition  to  variations  with  the  acoustic  frequency  and  the 
sediment  type,  sediment  attenuation  also  varies  with  depth  Into  the 
sediment.  Stoll Investigated  the  effects  of  depth  on  the  attenuation 
by  varying  the  overburden  pressure  In  his  theoretical  model.  He  found 
the  attenuation  decreases  with  Increasing  pressure  or  depth.  The 
variations  he  did  find,  however,  were  relatively  small  In  the  frequency 
range  between  10  and  1000  Hz,  decreasing  by  approximately  15X  over  a 
100  m  depth  change.  Laboratory  measurements  In  sand  sediments  confirm 
this  decrease,  although  the  attenuation  measurements  show  SS%  of  change 
for  the  same  depth  change.  For  the  high  porosity  sediments,  however, 
the  attenuation  measurements  reported  by  Hamilton  show  a  general 
increase  with  depth  In  the  first  400-600  m  of  sediment.  In  these 
sediments,  the  Increase  In  overburden  pressure  Is  accompanied  by  a 
decrease  In  the  porosity  as  the  depth  Increases.  The  change  In  the 
porosity  leaos  to  the  Increase  In  attenuation. 

There  are  two  viewpoints  on  sediment  attenuation,  that  differ 
on  how  attenuation  varies  with  frequency  and  depth.  Acoustic  analysis 
which  have  supported  the  Blot-Stoll  model  have  Ignored  the  possible 
Influence  of  the  depth  variations,  and  It  will  be  shown  later  how  this 
can  lead  to  erroneous  conclusions  about  the  frequency  dependence.  The 
point  of  view  taken  In  this  report  Is  guided  by  the  measurements.  The 
attenuation  Is  assumed  to  vary  linearly  with  the  acoustic  frequency  and 
to  vary  with  depth  Into  the  sediment. 

B.  Acoustic  Modeling  for  Shallow  Water 

The  estimation  procedure  developed  In  the  next  chapter  Is 
based  on  the  Interpretation  of  measured  propagation  loss  using  a 
propagation  model.  To  understand  the  procedure,  one  needs  to  know  how 
the  model  accounts  for  the  sediment  attenuation.  For  the  purpose  of 
estimating  environmental  parameters  from  acoustic  propagation  data, 
there  are  basically  two  ways  to  describe  the  acoustic  propagation:  ray 
theory  and  normal  mode  theory.  Both  theories  have  been  applied  to 
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numerical  modeling  of  the  propagation  In  shallow  water although  at 
frequencies  below  1000  Hz  the  normal  mode  theory  is  more  widely  used. 

Shallow  water  presents  difficulties  to  the  ray  theory 
approach.  In  the  vicinity  of  the  water-sediment  Interface,  the  basic 
ray  theory  assumptions  are  Invalid.^  In  shallow  water  the  entire  water 
column  Is  In  the  vicinity  of  this  Interface  for  the  low  frequencies. 
The  standard  procedure  around  this  limitation  Is  to  replace  the  sediment 
with  a  plane  wave  reflection  coefficient.  However,  at  low  frequencies  In 
shallow  water,  the  plane  wave  reflection  is  not  always  adequate,  and 


’•'•j 
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additional  corrections  are  needed  to  account  for  the  non-plane  wave 


nature  of  the  propagation. 
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Acoustic  propagation  In  shallow  water  can  be  expressed  by  the 
solution  to  the  separated  homogeneous  wave  equation  when  horizontal 
stratification  and  azimuthal  symmetry  are  assumed.  For  a  harmonic  time 
dependent  source,  the  acoustic  field  at  a  depth  z  and  range  r  from  the 
source  at  depth  z^  Is  expressed  as 


♦(r,z)  -  I  ♦„(*„)  ♦„(!)  H<’)(k„r) 


Where  H^^'(ICj^r)  Is  the  Hankel  function  of  the  first  kind,  Is  the 

nth  mode  or  eigenfunction  of  the  vertical  component  of  the  particle 
velocity,  and  k^  is  the  eigenvalue.  The  propagation  loss  can  be 
expressed  as 


ik  r,2l 


PL(r,z)  -  -10  109,o  fI  1?  *  "  1 


J*  1 
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where  the  asymptotic  form  of  H^^^(kj^r)  has  been  used.  The  1/r  term 
results  from  Hn)(|(^r)  and  represents  cylindrical  spreading  of  the 
modes.  Range  also  appears  In  the  exponents,  Ik^r.  Writing  the 
eigenvalue  k^  as 

k„  =  <n^^«n  •  n.4 
the  exponents  in  Eq.  II. 3  become 

attenuation  and  tc^r  Is  the  mode  phase.  An  Increase  In  range  Increases 
‘6j^r  for  each  n,  which  leads  to  a  general  Increase  in  PL(r,z).  However, 
small  changes  in  range  can  also  result  in  large  phase  differences 
between  the  modes.  This  can  result  In  constructive  or  destructive 
Interference  of  the  modes,  and  within  small  range  Intervals  PL(r,z)  can 
vary  as  much  as  15  dB.  These  fluctuations  In  PL(r,z)  can  easily  obscure 
the  general  trends  resulting  from  cylindrical  spreading  and  mode 
attenuation. 

Propagation  loss  measurements  are  often  made  using  explosive 
shots  and  are  evaluated  within  selected  frequency  bands.  The  frequency 
averaging,  represented  by  these  frequency  band  measurements,  eliminates 
most  of  the  fluctuations  In  the  propagation  loss  resulting  from  the  mode 
Interference.  These  measurements  can  usually  therefore  be  characterized 
by  cylindrical  spreading  and  attenuation. 

The  mode  Interference  can  easily  be  removed  from  the 
calculations  (Eq.  II. 3)  by  expressing  the  propagation  loss  as  an 
Incoherent  sum  of  the  modes 


PL(r.z) 


-10  log^Q 


II. 5 
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Equation  II. 5  defines  propagation  loss  as  a  monotonically  Increasing 
function  of  range  which  can  be  characterized  by  cylindrical  spreading 
and  mode  attenuation. 


The  calculation  of  acoustic  propagation  now  Involves  the 

determination  of  k.  and  0..  If  shear  waves  are  not  considered,  the 

n  n 

depth  component  of  the  wave  equation  takes  the  form 


II. 6 


where  k(z)  Is  the  propagation  wave  number.  The  solutions  to  Eq.  II. 6 
must  also  satisfy  the  boundary  conditions 


♦„  (0)  =  0 


9Z 


3Z 


II. 7 


p(z)i|»^(2) 


=  p(z)*^(z) 


where  the  H^'s  designate  the  horizontal  boundaries.  These  conditions 
result  from  a  pressure  release  surface  condition,  and  continuity  of 
pressure  and  normal  component  of  particle  velocity  at  each  horizontal 
boundary  Identified  by  a  change  In  density.  The  Sonmerfeld  radiation 
condition  Is  also  Imposed  on  the  solution  to  allow  only  outward 
propagation  at  Infinite  depth. 


Equation  II. 6  will  take  into  account  attenuation  in  the  medium 
when  k(z)  is  defined  as  a  complex  function 

k{z)  =  +  ia(z)  ,  II. 8 

where  a(z)  is  the  attenuation,  c(z)  is  the  sound  speed,  and  u  is  the 
acoustic  angular  frequency.  In  this  case  both  0^(z)  and  k^  will  also  be 
complex.  The  mode  attenuation  coefficient  6^  is  defined  as  the  imaginary 
part  of  k^,  as  in  Eq.  I I. 4. 

Following  Schiff^^  and  Koch  et  al.,^^  a  perturbation  solution 
for  6^  can  be  developed  in  terms  of  and  <^.  This  approach  is 
based  on  the  assumption  that  a(z)«u/c(z).  The  attenuation  is  introduced 
as  a  slight  perturbation  to  the  environment  with  the  effect  of  producing 
a  slight  perturbation  to  k^  and  Equation  II. 6  is  rewritten  by 
substituting  Eq.  II. 8  for  k(z),  Eq.  II. 4  for  k^,  (k^+Ak^)  for  and 
((/y,+A^)  for  0^  .  This  new  Eq.  II. 6  is  expanded  and  all  terms  with  two 
or  more  perturbation  parameters  are  dropped.  This  final  differential 
equation  is  linear  with  respect  to  the  perturbation  parameters.  By 
integrating  this  equation,  6^  can  be  expressed  as 


''  II. 9 

^n“  dz  , 

where  1^  is  the  normalization  factor  for  the  nth  mode.  The  perturbed 
solutions  for  k  and  <|i„(z)  are  the  same  as  those  calculated  without 

n  rip^  pr 

considering  attenuation.  The  mode  attenuation  coefficients  6^  are 

therefore  expressed  in  terms  of  the  zero  attenuation  solutions. 
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Sediment  attenuations  can  be  handled  by  solving  Eq.  II. 6  without 
attenuation,  calculating  6^,  by  Eq.  II. 9,  and  defining  a  new  eigenvalue 
by  Eq.  II. 4. 

C.  Effects  of  Depth  Variations 

The  effects  of  the  fact  that  the  attenuation  varies  with  depth 
can  be  examined  by  looking  at  Eq.  II. 9.  When  the  sediment  attenuation 
is  assumed  to  be  constant  with  depth,  Eq.  I I. 9  simplifies  to 

*^n  ~  “sed^nl  ^  “sub^n2  * 


_  Psed^  f 

•'n^n  •^ed 

^sub*^  fit  19  1 

iubi  ■=<"> 

where  the  subscripts  "sed"  and  "sub"  denote  properties  of  the  bottom 
sediment  and  the  substrate  underlying  tbe  sediment,  respectively. 
Except  for  thin  sediments,  the  contribution  from  the  substrate  y^2 
be  ignored  for  the  low  order  modes.  Experimental  measurements  of  mode 
attenuation  coefficients  S  can  be  used  to  estimate  the  attenuation  in 
the  sediments^®"^®*^®*^^  by 


♦„(Z) 


c(z)‘^  dz 


II. 10 


a'(f)  =  S^/y 


nl 


II. 11 


Equation  II. 11  can  result  in  attenuation  estimates  which  are 
not  linear  with  frequency.  This  can  be  illustrated  for  a  hypothetical 
environment  by  calculating  the  mode  attenuation  coefficient  according  to 
Eq.  II. 9  and  then  estimating  the  sediment  attenuation  by  Eq.  II. 10.  Two 
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attenuation  profiles,  typical  of  clays,  were  used  as  examples 

(Fig.  II. 1).  Both  have  a  value  of  0.03  dB/m-kHz  in  the  first  2  m,  and 
yet  they  yield  significantly  different  power  laws  for  6^  in  the  50- 
800  Hz  region.  The  coefficients  6^  were  computed  using  Eq.  II. 9  at 
frequencies  of  50,  100,  200,  400,  and  800  Hz  with  a  linear  relation  of 
attenuation  to  frequency.  The  required  's  and  k.'s  were  calculated 

pO  OQ  ^ 

using  the  computer  program  NEMESIS.  *  The  6^^'$  were  treated  as  if 
they  were  measured  data  and  converted  to  equivalent  constant-depth 

attenuation  a'(f)  by  Eq.  II. 11.  The  results  are  presented  by  the  points 
in  Fig.  II. 2.  Also,  the  actual  attenuation  at  the  sediment  surface  is 

shown  by  the  solid  line.  The  least  square  error  curves  (indicated  by 

the  dashed  lines  in  Fig.  I I. 2)  for  these  points  were 

=  29.5  f^*^^dB/km  ,  11.12 

with  a  mean  square  error  of  0.09  (dB/km)  for  profile  1,  and 

*2  =  33.0  f°’^^  dB/km  ,  11.13 

with  a  mean  square  error  of  0.03  (dB/km)  for  profile  2.  The 
frequencies  are  in  kilohertz  in  Eqs.  11.12  and  11.13. 

The  frequency  dependence  of  a',  determined  by  Eq.  1 1. 11,  can 
be  strongly  affected  by  the  depth  dependence  of  the  actual  attenuation. 
For  an  attenuation  profile  which  increases  with  depth,  «'  at  the  lower 
frequencies  is  greater  than  the  attenuation  at  the  sediment  surface. 
This  would  lead  to  n  in  Eq.  II. 1  being  less  than  1.  When  the 
attenuation  profile  decreases  with  depth,  a'  at  the  lower  frequencies  is 
less  than  the  attenuation  at  the  sediment  surface.  This  would  lead  to  n 
being  greater  than  1. 
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FIGURE  1L2 

SEDIMENT  ATTENUATIONS  a'(f)  CALCULATED  USING  Eq.  0.11 
FOR  THE  ATTENUATION  PROFILES  GIVEN  IN  Fig.  0.1 
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Ingenito  found  n=1.75  while  Rubano  found  n=l.l,  which  « 

would  suggest  attenuation  profiles  which  decrease  with  depth.  Ingenito  | 

worked  In  a  region  with  a  sandy  bottom.  It  has  already  been  mentioned  j 

that  attenuation  Is  expected  to  decrease  with  depth  for  a  sandy  bottom. 


The  environmental  description  discussed  In  Ref.  17  and 
references  therein  were  used  along  with  the  attenuation  In  Fig.  II. 3  to 
compute  the  mode  attenuation  coefficients  at  400  and  750  Hz.  A  linear 
relationship  between  attenuation  and  frequency  was  also  assumed. 
Equation  II. 11  was  used  to  compute  a*  for  each  frequency,  which  was  then 
compared  to  the  values  reported  by  Ingenito  (Fig.  I 1.4).  The  calculated 
values  result  In  n»1.22.  Slight  changes  could  be  made  In  the 
attenuation  near  the  surface  to  Increase  the  mode  attenuation  at  750  Hz 
while  leaving  400  Hz  virtually  unchanged,  thus  Increasing  the  apparent 
power  law.  However,  the  calculated  values  do  fall  within  the  error  bars 
reported  with  the  measurements  showing  that  an  attenuation  profile, 
linearly  dependent  on  frequency,  can  explain  the  measurements. 

The  frequency  dependence  reported  by  Rubano  closely  matches 
the  results  In  Fig.  II. 2  for  profile  1.  Using  the  dashed  curve  In 
Fig.  II. 1  to  define  an  attenuation  profile  and  the  environmental 
description  In  Ref.  16,  a's  were  computed  at  50,  100,  and  400  Hz. 
Figure  II. 5  presents  a  comparison  between  these  calculations  and  the 
results  reported  by  Rubano.  The  largest  difference  was  0.35  dB/km  at 
400  Hz.  Here  again,  an  attenuation  profile,  linearly  dependent  on 
frequency,  can  explain  the  measured  results. 

Propagation  loss  measurements,  such  as  those  of  Rubano  and 
Ingenito,  can  provide  estimates  of  attenuation  profiles  In  the  sediment 
which  are  linearly  dependent  on  frequency.  In  the  next  chapter  a 
technique  for  estimating  these  profiles  Is  developed  based  on 
propagation  loss  measurements  collected  at  numerous  frequencies. 
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FIGURE  IL6 

COMPARISON  OF  CALCULATED  DEPTH-CONSTANT 
ATTENUATIONS  FOR  THE  MODIFIED  VERSION  OF 
PROFILE  1,  Fi9.1L1,  AND  RESULTS  FROM  R«f.  16 
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CHAPTER  III 

ESTIMATING  THE  AHENUATION  PROFILE 

The  Importance  of  the  depth  variations  In  the  attenuation  was 
brought  out  In  the  discussions  presented  In  the  preceding  chapter.  When 
acoustic  propagation  data  are  analyzed  under  the  assumption  that 
attenuation  Is  Invariant  with  depth,  one  can  misinterpret  the 
attenuation  dependence  on  frequency  or  overemphasize  the  Importance  of 
other  loss  mechanisms.  To  avoid  these  pitfalls,  the  attenuation  needs 
to  be  described  and  analyzed,  not  only  as  a  function  of  frequency,  but 
also  as  a  function  of  depth.  A  method  or  procedure  Is  needed  which  will 
estimate  the  attenuation  depth  profile  within  a  measurement  region. 

A  method  can  be  devised  for  estimating  the  attenuation  profile 
which  will  take  advantage  of  variations  In  the  shape  of  the  low  order 
normal  modes  as  the  frequency  varies.  These  variations,  combined  with 
the  frequency  dependence  of  the  sediment  attenuation,  result  In  the 
observed  frequency  dependence  In  the  measured  propagation  loss. 
Estimates  can  be  obtained  without  requiring  direct  In  situ  measurements 
of  the  sediment  properties.  The  only  measurements  required  are  those 
generally  taken  during  acoustic  propagation  experiments. 

A.  Input  Parameters 

The  estimation  technique  requires  propagation  loss 
measurements  and  an  environmental  description.  The  environmental 
description  Includes  the  sound  speed  profile  In  both  the  water  and 
sediment,  the  density  of  both  the  water  and  the  sediment,  and  an 
estimate  of  the  thickness  of  the  sediment.  Propagation  loss 
measurements  collected  at  multiple  frequencies  are  affected  by  the 
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acoustic  attenuation  In  the  sediment  and  can  be  used  to  estimate  the 
attenuation  profile. 


Sediment  attenuation  values  are  assigned  to  various  depths 
based  on  characteristics  of  the  normal  modes  to  be  discussed  below.  The 
calculations  of  the  normal  modes  require  Input  parameters  defining  the 
acoustic  environment.  The  sound  speed  profiles  In  both  the  water  and 
the  sediment  layers  are  the  most  Important  of  these.  The  sound  speed 
profiles  enter  Into  the  calculation  of  the  modes  through  Eq.  11.6,  and 
thus  determine  how  much  energy  will  penetrate  Into  the  bottom  sediments. 
The  sound  speed  In  the  sediment  also  enters  Into  the  determination  of 
the  mode  attenuation  term  through  Eq.  II. 9.  Changes  to  the  sound  speed 
profile  will  result  In  changes  In  the  calculated  modes  and  In  the 
calculated  mode  attenuation  terms. 

Propagation  loss  measured  as  a  function  of  range  Is 
significantly  Influenced  by  the  attenuation  In  the  sediment. 
Information  about  the  attenuation  Is  contained  within  these 
measurements.  The  propagation  loss  has  to  be  measured  at  several 
frequencies  In  order  to  estimate  the  attenuation  profile.  There  Is  a 
sediment  layer  defined  for  each  frequency,  and  therefore  the  depth 
resolution  of  the  estimates  Is  a  function  of  the  number  of  frequencies 
for  which  there  are  measurements. 

B.  Frequency  Characteristics  of  the  Normal  Modes 

There  Is  one  characteristic  of  the  normal  modes  that  Is 
Important  to  the  procedure  for  estimating  the  attenuation  profile  In 
shallow  water  sediments.  This  characteristic  will  be  referred  to  as  the 
penetration  depth  of  the  mode.^  This  depth  defines  how  much  of  the 
sediment  affects  the  mode.  This  depth  has  also  been  called  the  hidden 
depth,  meaning  that  the  sediment  below  this  depth  Is  hidden  from  the 
mode.  For  any  specific  sediment,  this  penetration  depth  Is  a  function 
of  frequency.  As  the  frequency  decreases,  the  penetration  depth 
Increases  Into  the  sediment.  This  frequency  dependence  can  be 
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demonstrated  by  looking  at  the  results  obtained  In  a  shallow  water 

environment  consisting  of  two  Isovelocity  layers:  c^  Is  the  sound  speed 

In  the  water  layer  of  thickness  H,  and  C2  Is  the  sound  speed  of  the 

sediment  layer  of  semi-infinite  thickness.  This  was  the  description 
30 

used  by  Pekerls,  who  noted  the  frequency  dependence  of  the  penetration 
depth.  For  this  environment,  the  mode  function  in  the  sediment  layers 
Is  described  by  decaying  exponentials. 


^„(z)  «  D  exp  (kI  -  zj. 


III.l 


where  z  Is  depth,  k  «  ^/C2>  0  Is  a  constant  satisfying  the  boundary 
conditions,  is  the  eigenvalue,  and  j  *  the  angular  frequency, 
approaches  lj/c^  at  the  higher  frequencies  and  approaches  lj/c2  at  the 
lower  frequencies,  Is  monotonic  In  ij  between  these  two  limits.  As 
the  frequency  Increases,  the  decay  constant  In  Eq.  III.l, 

Increases,  and  the  mode  function  decreases  more  rapidly  with  depth. 

In  Chapter  II,  the  mode  attenuation  term  was  found  to  be 


’"'Vn  till 


This  equation  shows  that  the  Influence  of  the  attenuation  within  a 
sediment  layer  Is  dependent  upon  the  Integral  across  that  layer,  of  the 
eigenfunction  squared.  To  understand  the  interaction  between  the  mode 
and  the  sediment  attenuation,  consider  a  hypothetical  attenuation 
profile  defined  by 

a(z)  *  a  +  g(z) 

with  the  mode  function  In  the  sedlnwnt  given  by  Eq.  III.l.  The  partial 
contribution  to  6^  from  the  sediment  at  depths  less  than  z  Is  given  by 


d„(z)  -  [ 

1 

g  +  Eia^®  J^n 

> 

III. 2  ; 
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Mhere  At  a  fixed  z,  as  the  frequency  Increases 
(Increasing  T  ),  the  second  term  In  Eq.  I I I. 2  decreases,  and  therefore 
the  relative  contribution  Increases.  The  shallower  sediment  layers  will 
contribute  more  to  the  total  mode  attenuation  at  the  higher  frequencies 
while  the  deeper  sediment  layers  will  only  contribute  to  the  lower 
frequency  attenuations. 

When  the  sound  speed  In  the  sediment  changes  with  depth, 
Eq.  III.l  still  approximates  the  low  order  mode  functions  within  the 
sediment.  The  Influence  of  the  deeper  sediments  upon  mode  attenuation 
Is  therefore  expected  to  Increase  with  decreasing  frequency.  This  can 
be  demonstrated  by  looking  at  the  first  mode  for  two  frequencies,  50  and 
400  Hz,  for  more  realistic  environments.  The  first  environment  Is  430  m 
of  Isovelocity  water  overlying  50  m  of  a  clay  sediment  (Fig.  III.l). 
The  sediment  sound  speed  Is  defined  by  an  Interface  sound  speed  ratio  of 
1.0  at  the  water-sediment  boundary,  and  by  a  sound  speed  gradient  of 
1,0  sec“^.  This  Is  a  nominal  description  for  a  clay  sediment. 
Figure  III.l  presents  mode  1  at  both  50  and  400  Hz  for  the  environment 
described  above.  At  50  Hz  the  mode  extends  down  through  the  sediment 
layer,  while  at  400  Hz  an  appreciable  amplitude  extends  down  only  to 
10  m.  The  Influence  of  attenuation  at  depths  below  10  m  Is  expected  to 
be  greater  at  50  Hz  than  at  400  Hz.  Figure  1 1 1. 2  presents  a  more 
detailed  look  at  the  modes  within  the  sediment.  The  dotted  curves  are 
the  mode  functions  within  the  sediment  normalized  to  a  value  of  one  at 
the  sediment  surface.  The  mode  function  at  50  Hz  decreases  more  slowly 
than  the  one  at  400  Hz,  as  suggested  by  Eq.  III.l. 

As  was  stated  above,  the  mode  attenuation  coefficient  Involves 
the  Integral  of  the  mode  function  squared.  The  solid  curves  In 
Fig.  III. 2  present  the  accumulative  value  for  this  Integral  as  a 
function  of  depth  Into  the  sediment.  These  curves  have  been  normalized 
so  that  the  total  Integral  has  a  value  of  one,  and  are  related  to 
Eq.  III. 2.  At  400  Hz,  98%  of  this  Integral  occurs  above  10  m.  At 
50  Hz,  however,  only  55%  of  the  Integral  has  occurred  within  this  10  m. 
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ENVIRONMENTAL  DESCRIPTION  AND  MODE  1  AT  50  Hi 
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The  sediment  attenuation  in  the  first  10  m  almost  totally 
determines  the  mode  attenuation  coefficient  at  400  Hz.  Sediment 
attenuation  values  down  to  40  m  are  important  in  che  determination  of 
the  coefficient  at  50  Hz,  however. 

A  second  environment  to  consider  is  90  m  of  water  overlying 
30  m  of  a  sand  sediment.  The  sound  speed  profile  in  the  water  is 
downward  refracting,  and  the  sound  speed  profile  in  the  sand  follows 
descriptions  given  by  Hamilton.^  Figure  III. 3  presents  the  normalized 
mode  function  and  mode  function  integral  in  the  sediment  for  50,  400, 
and  1600  Hz.  From  these  figures,  it  is  seen  that  the  modes  in  this  case 
have  depth  dependencies  similar  to  those  in  the  previous  case.  The 
higher  frequencies  only  penetrate  to  shallowest  depths.  The  lower 
frequencies  penetrate  the  deeper  sediments,  although  the  depths  of 
penetration  are  less  than  those  seen  in  Fig.  III. 2.  The  mode 
attenuation  for  the  highest  frequencies  are  determined  by  just  the 
surface  sediments,  and  as  the  frequencies  decrease  more  of  the  sediment 
enters  into  this  determination. 


Estimating  the  Attenuation  Profile 


The  relationship  between  the  frequency  and  the  penetration 
depth  presented  above  suggests  a  method  for  estimating  the  sediment 
attenuation  profile.  By  examining  the  frequency  dependence  of  the 
propagation  loss  in  terms  of  the  penetration  depths,  one  can  estimate  an 
attenuation  profile. 


An  estimation  procedure  was  developed  around  the  relationship 
between  frequency  and  penetration  depth.  It  uses  propagation  loss 
measured  at  multiple  frequencies  and  the  linear  frequency  dependence  of 
the  attenuation  in  the  sediment.  The  sediment  can  be  divided  into 
layers  according  to  the  penetration  depths  for  the  frequencies  available 
in  the  measurements.  Estimates  of  the  sediment  attenuation  within  each 
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layer  are  determined  from  the  propagation  loss  measurements  at  the 
frequency  corresponding  to  that  layer. 

The  procedure  is  divided  into  two  parts,  as  illustrated  in 
Fig.  1 1 1. 4.  The  first  step  is  to  define  layer  depths  and  the  second  is 
to  determine  the  attenuation  values  in  each  layer.  Propagation  beyond 
a  few  water  depths  involves  only  the  lowest  order  modes.  Therefore,  for 
each  frequency  one  sediment  layer  will  be  defined  by  considering  only 
mode  1.  The  integral  of  the  mode  function  squared  is  computed  as  a 
function  of  depth,  as  presented  in  Figs.  III. 2  and  III. 3.  A  layer 
depth  can  now  be  defined  as  the  depth  at  which  this  integral  exceeds  a 
specified  fraction  of  the  integral  computed  across  the  entire  sediment 
column.  For  example,  using  80%  of  the  total  integral,  layer  depths  of 
4.7  m  and  18  m  are  obtained  for  the  modes  presented  in  Fig.  III. 2,  and 
depths  of  0.5  m,  1.2  m,  and  8.6  m  for  Fig.  III. 3.  If  there  are  data  at 
n  frequencies  available,  there  will  be  n  sediment  layers  defined.  The 
last  layer,  instead  of  extending  down  to  the  80%  depth  of  the  lowest 
frequency,  would  extend  down  to  the  bottom  of  the  sediment  column. 

With  the  layer  depths  defined,  the  attenuation  values  can  be 
estimated  by  comparing  the  measured  propagation  losses  to  the  calculated 
losses.  The  estimates  of  sediment  attenuation  a  are  initially  assumed 
to  be  the  same  in  all  the  sediment  layers.  Then,  beginning  with  the 
highest  frequency  of  interest  for  which  there  are  measurements,  800  Hz 
in  Fig.  III. 4,  the  sediment  attenuation  a  is  varied  to  minimize  the 
difference  between  the  measured  and  calculated  results.  The  calculated 
propagation  losses  are  computed  using  the  incoherent  mode  sum  given  in 
Eq.  II. 5.  The  value  thus  obtained,  is  assigned  to  the  initial 
sediment  layer  and  varies  linearly  with  frequency.  The  model/data 
comparisons  are  now  conducted  for  the  next  lower  frequency  of  interest. 
The  attenuation  in  the  sediment  below  the  first  layer  is  again  assumed 
to  be  constant  between  these  deeper  layers  and  is  varied  to  minimize  the 
propagation  loss  differences  at  this  frequency.  The  attenuation  value 
obtained,  "fs  assigned  to  the  second  layer,  again  varying  linearly 
with  frequency.  The  procedure  is  now  repeated  for  the  third 
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ATTENUATION  PROFILE  GENERATED  BY  FREQUENCY  CONSIDERATIONS 
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highest  frequency,  and  is  continued  until  all  available  frequencies  have 
been  considered.  The  final  result  is  an  attenuation  profile,  described 
by  n  layers  of  constant  attenuation  and  a  linear  frequency  dependence. 

Calculated  propagation  losses  versus  range  were  used  to 
simulate  measured  values  based  on  a  hypothetical  environment.  In  this 
way,  the  actual  attenuation  profile  was  known,  and  a  comparison  between 
actual  and  estimated  profiles  can  be  performed.  The  environment 
presented  in  Fig.  III.l  along  with  the  two  attenuation  profiles 
represented  by  the  solid  curves  in  Fig.  1 1 1. 5  were  used  to  compute 
propagation  losses.  Propagation  losses  were  calculated  as  a  function  of 
range  at  frequencies  of  50,  100,  200,  400,  and  800  Hz  using  the  computer 

OO  pQ 

normal  mode  model  NEMESIS.  *  These  calculated  values  of  propagation 
loss  were  then  used  as  if  they  were  data  in  the  estimation  procedure 
outlined  in  Fig.  III. 4.  The  estimates  were  based  on  the  t//^(z)'s  and 
K  's  obtained  for  the  propagation  loss  calculations. 

The  estimated  attenuation  profiles  are  presented  in  Fig.  III. 5 
as  dashed  lines.  For  these  estimates,  the  sediment  layer  depths  were 
defined  by  80%  of  the  total  mode  square  integral.  In  both  cases  the 
five-layer  profile  approximates  the  input  profile,  with  the  major 
difference  occurring  in  the  first  layer.  The  first  layer  as 
attenuations  are  biased  slightly  towards  the  deeper  layer  values. 
Figure  III. 6  presents  the  propagation  loss  comparisons  for  the  actual 
and  the  estimated  profiles.  The  maximum  error  between  any  two 
corresponding  curves  is  less  than  1  dB  and  the  maximum  mean  square  error 
for  any  one  frequency  was  0.2  dB. 

The  comparison  just  presented  was  an  ideal  case.  The 
environment  was  known  exactly  and  the  propagation  loss  inputs  were  exact 
and  not  corrupted  by  noise.  For  actual  measurements,  the  environmental 
description  will  generally  be  composed  of  several  measured  parameters, 
such  as  the  water  sound  speed  profile,  while  the  remaining  parameters 
must  be  inferred  from  indirect  methods.  The  sediment  sound  speeds  and 
densities  are  inferred  from  the  sediment  type  and  subbottom  profiling. 
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CHAPTER  IV 

AHENUATION  PROFILES  IN  EXAMPLE  SHALLOW  WATER  REGIONS 


Chapter  III  presented  a  procedure  for  estimating  the 
attenuation  profile  In  shallow  water  regions  which  Is  based  on  acoustic 
propagation  measurements.  The  procedure  will  be  applied  to  measured 
propagation  loss  data  collected  In  three  regions.  Data  collected  In 
each  region  provided  the  Information  necessary  to  obtain  the  estimates. 
Acoustic  propagation  measurements  were  obtained  using  explosive  sources, 
that  generated  broadband  signals.  This  provided  propagation  loss 
measurements  at  multiple  frequencies,  as  required.  There  was  also 
sufficient  Information  collected  about  the  environments  to  allow  the 
sound  speed  profiles  In  both  the  water  and  sediment  layers  to  be 
described.  The  normal  modes  could  therefore  be  calculated  for  each 
region. 

A.  Environmental  Descriptions 

Figures  IV.1«IV.3  present  the  measured  propagation  loss  and 
the  bathymetry  along  the  source  tracks  for  three  measurement  sites. 
These  three  sites  provided  a  variety  of  environmental  conditions:  three 
different  water  depths,  two  sediment  types,  and  two  water  sound  speed 
profiles. 


Site  1,  Fig.  IV. 1,  lies  within  a  sedimentary  basin.  The 
sediment  In  this  region  of  the  basin  has  been  classified  as  clay  based 
on  seismic  refraction  profiling  and  sediment  coring  samples.  The  sound 
speed  profile  was  defined  by  a  water-sediment  sound  speed  ratio  of  1,  a 
gradient  of  1  sec~^,  and  a  thickness  of  50  m.  The  water  depth  at  this 
site  was  430  m  and  the  sound  speed  profile  was  described  as  Isovelocity 
(Fig.  III.l). 
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BATHYMETRY  AND  MEASURED  PROPAGATION  LOSS  FOR  SITE  1 
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BATHYMETRY  AND  MEASURED  PROPAGATION  LOSS  FOR  SITE  2 
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Site  2,  in  Fig.  IV. 2,  was  located  250  km  from  site  1  on  a 
bathymetric  high  outside  the  sedimentary  basin  at  site  1.  The  water 
depth  was  250  m  and  again  the  sound  speed  profile  was  isovelocity. 
Based  on  seismic  profiling  conducted  in  this  area,  the  sediment  at  this 
site  was  assumed  to  be  50  m  of  clay  with  the  same  description  as  site  1. 
The  attenuation  in  the  sediment  at  this  site,  however,  was  not  expected 
to  be  the  same  as  at  site  1. 

Site  3,  Fig.  IV. 3,  was  located  on  the  continental  shelf  in 
120  m  of  water.  The  water  sound  speed  profile  at  this  site  was  downward 
refracting,  sound  speed  decreasing  with  depth.  Sediment  coring  samples 
and  seismic  refraction  profiling  indicated  the  sediment  type  was  silty 
sand.  Within  the  first  meter  of  sediment,  the  sound  speed  increased  to 
1670  m.  The  sound  speed  gradient  decreased  with  depth  throughout  the 
30  m  of  sediment. 

B.  Measured  Propagation  Loss 

Figures  IV. 1-3  also  present  the  propagation  losses  versus 
range  for  each  site.  The  sources  were  detonated  at  a  depth  of  91  m  for 
all  three  sites.  The  receiver  depths  were  91  m  at  sites  1  and  2,  and 
75  m  at  site  3. 

The  propagation  losses  to  all  three  sites  were  relatively  low. 
The  losses  at  most  frequencies  did  not  exceed  95  dB  at  sites  1  and  2,  or 
90  dB  at  site  3,  even  out  to  ranges  of  140  km. 

At  site  1  the  measured  propagation  losses  in  the  range 
interval  0-80  km  was  used  to  estimate  the  attenuation  profile.  The 
selection  of  this  range  interval  was  based  on  the  bathymetry.  Beyond 
80  km,  there  was  a  decrease  in  water  depth.  To  avoid  the  effects  of 
change  in  water  depth,  the  analysis  was  confined  to  the  first  80  km. 


At  site  2,  Fig.  IV. 2,  the  propagation  losses  show  a  change  in 
the  range  dependence  at  a  range  of  approximately  55  km.  This  change 


correlates  with  a  known  change  in  the  sediment,  and  therefore  only 
ranges  less  than  55  km  were  used  in  the  analysis.  Within  this  range 
interval  the  propagation  exhibited  little  dependence  on  frequency. 

At  site  3  there  is  no  obvious  change  in  the  range  dependence 
of  propagation  as  seen  at  site  2.  As  at  site  1,  there  was  a  change  in 
water  depth.  The  depth  changed  from  120  m  at  site  3  to  160  m  at  the  end 
of  the  source  track.  The  analysis  was  limited  to  the  first  100  km  to 
avoid  possible  effects  due  to  the  changing  water  depth. 

C.  Attenuation  Profiles 

The  attenuation  profiles  were  estimated  for  each  of  the  three 
sites,  based  on  the  data  presented  in  the  last  two  sections.  The 
estimated  attenuation  values  were  all  within  the  ranges  expected  for  the 
sediment  type  found  at  the  measurement  sites.  The  depth  dependence  of 
each  profile  also  agreed  with  previously  reported  measurements. 

Figure  II.l  presented  the  attenuation  profiles  estimated  for 
sites  1  (curve  1)  and  2  (curve  2).  The  decrease  in  attenuation  at  site  1 
was  needed  to  account  for  the  lower  losses  at  the  lower  frequencies. 
The  increase  in  attenuation  with  depth  at  site  2  was  needed  to  account 
for  the  near  frequency  independent  nature  of  the  propagation  loss.  The 
differences  between  these  two  profiles  indicate  the  variability  in  the 
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attenuation  values  for  clay  type  sediments.  Data  compiled  by 
Hamilton^*^  had  indicated  that  attenuation  could  either  increase  or 
decrease  with  increasing  depth. 

A  comparison  between  the  data  and  the  final  calculations  at 
site  1  showed  that  the  calculated  losses  generally  were  within  3  dB  of 
the  measurements.  The  comparisons  at  50  Hz  and  400  Hz  are  shown  in 
Fig.  IV. 4.  Comparisons  at  the  other  frequencies  presented  similar 
results.  The  rms  differences,  calculated  at  each  frequency,  varied  from 
2.3  dB  at  400  Hz  to  1.6  dB  at  100  Hz. 


Comparisons  of  the  final  calculations  and  the  measurements  at 
site  2  were  similar  to  those  found  at  site  1.  The  calculations  for  this 
site  were  also  generally  within  3  dB  of  the  measurements  with  the  rms 
differences  varying  from  1.7  dB  at  400  Hz  to  1.4  dB  at  50  Hz.  Fig¬ 
ure  IV. 5  presents  the  comparisons  at  50  Hz  and  400  Hz.  Comparisons  at 
other  frequencies  were  similar,  with  the  curves  lying  between  those  at 
50  and  400  Hz. 

For  both  sites  1  and  2,  where  the  sediments  were  clay,  the 
attenuation  values  were  In  the  lower  range  of  values  expected  for  clays. 
These  estimates  varied  from  0.01  up  to  0.08  dB/m-kHz.  Although  the 
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attenuation  Is  generally  expected  to  Increase  with  depth,  the  fact  that 
at  site  1  the  profile  decreases  with  depth  Is  not  a  point  of  contention. 
The  expected  Increase  Is  based  on  a  compilation  of  reported  results 
within  the  first  400-500  m.  The  slight  decrease  In  attenuation  within 
the  first  50  m  of  sediment  Is  not  a  significant  deviation  from  a  trend 
seen  In  such  a  diversified  set  of  measurements. 

The  estimated  profile  obtained  for  the  sand  sediment  was  also 

g 

In  agreement  with  the  nominal  profile  (Fig.  IV. 6).  The  attenuation 
varied  from  0.28  to  0.6  dB/m-kHz,  which  Is  In  the  middle  of  the  range 
for  sands.  The  attenuation  decreased  with  Increasing  depth.  This  was 
similar  to  the  typical  sand  profile  presented  earlier  In  Fig.  II. 3. 

The  comparisons  between  calculated  and  measured  propagation 
losses  at  site  3  were  similar  to  those  found  at  sites  1  and  2, 
Fig.  IV. 7.  The  calculations  were  generally  within  3  dB  of  the 
measurements  and  the  rms  differences  varied  from  2.8  dB  at  400  Hz  to 
3.5  dB  at  50  Hz.  The  larger  rms  errors  at  this  site  were  probably  due 
to  an  Increase  In  the  point  to  point  variation  of  the  measurements. 

All  three  profiles  estimated  by  the  procedure  in  Chapter  III 
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were  In  general  agreement  with  measurements  reported  by  Hamilton.  The 
rms  differences  In  the  calculated  propagation  losses  were  significantly 
larger  than  those  found  In  the  test  case  at  the  end  of 
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Chapter  III,  2  dB  compared  to  0.2  dB.  There  are  several  possible 
sources  of  these  errors.  The  estimated  profiles  were  based  on  sediment 
sound  speeds  derived  strictly  from  the  sediment  tyf>e.  The  measured 
propagation  loss  did  not  strictly  conform  to  the  Incoherent  summation  of 
the  modes  assumed  by  the  calculations.  The  data  displayed  point  to 
point  variations  resulting  from  possible  experimental  errors,  variations 
1n  the  noise,  and  variations  In  the  actual  propagation  losses.  Each  of 
these  deviations  lead  to  the  Increased  differences  found  between  the 
calculated  and  measured  propagation  losses. 

D.  Conclusions  and  Recommendations 

The  examples  above  show  that  the  attenuation  profiles  In 
sediments  of  the  bottoms  of  shallow  ocean  regions  can  be  estimated  using 
propagation  loss  measurements.  The  estimation  technique  In  Chapter  III 
Is  based  on  normal  mode  theory  and  the  fact  that  propagation  beyond  a 
few  water  depths  Involves  only  the  lower  order  modes.  The  depth  of 
penetration  of  these  modes  are  frequency  dependent.  The  low  frequencies 
penetrate  deeper  than  the  higher  frequencies.  The  frequency-depth  of 
penetration  relationship  and  the  frequency  dependence  of  the  measured 
propagation  loss  are  used  to  obtain  the  attenuation  profile  estimates. 

The  estimated  attenuation  profiles  discussed  In  this  chapter 
(as  shown  In  Figs.  II, 1  and  III.4),  agree  with  the  characteristics  of 
attenuation  profiles  measured  In  situ  and  In  the  laboratories,^ 
discussed  In  Chapter  II.  The  agreements  are  between  the  estimates  and 
measurements  made  In  similar  sediment  types.  The  estimates  In  clays  are 
In  the  range  0.01-0.1  dB/m-kHz,  typical  of  clays.  In  sand  the  estimates 
are  In  the  range  0.2-0. 9  and  decrease  with  depth,  which  Is  typical  of 
attenuation  In  sand. 

The  agreements  found  between  these  attenuation  profiles  and 
laboratory  and  In  situ  measurements  Is  encouraging.  However,  there  are 
several  recommendations  that  can  Improve  the  estimation  procedure.  In 
the  examples  presented  above,  only  a  limited  portion  of  the  data  sets 


were  used  in  the  estimation  of  attenuation.  Estimation  beyond  some 
range  was  eliminated  because  of  a  significant  change  in  the  environment. 
It  should  be  possible  to  expand  and  modify  the  present  technique  to  use 
the  longer  range  data  in  obtaining  a  range  dependent  attenuation 
estimate.  Besides  extending  the  usefulness  of  the  propagation  data, 
this  would  provide  a  more  complete  description  of  a  measurement  region. 

Another  step  that  is  needed  is  an  error  analysis  of  the 
attenuation  estimates.  The  calculated  incoherent  propagation  loss 
curves  are  smooth  monotonic  functions  of  range.  The  measured  losses, 
however,  contain  variations  about  a  general  trend  of  increasi’ig  loss 
with  increasing  range.  These  variations  result  from  interference 
between  modes,  inhomogeneities  in  the  environment,  and  experimental 
error.  These  variations  in  the  measured  losses  introduce  uncertainties 
in  the  attenuation  estimates.  More  sophisticated  calculations  of 
propagation  loss  could  account  for  the  interference  between  modes.  The 
variations  due  to  the  inhomogeneities  and  measurement  errors,  however, 
will  still  lead  to  uncertainties  in  the  estimates.  At  the  moment,  the 
estimation  technique  provides  an  attenuation  profile,  but  does  not 
indicate  what  the  confidence  limits  or  the  error  bounds  are. 

The  difficulty  in  setting  the  confidence  limits  is  deciding 
what  those  limits  mean  and  how  they  should  be  determined.  There  will  be 
errors  associated  with  attenuation  values  assigned  to  each  sediment. 
There  will  also  be  errors  associated  with  the  fact  that  the  sediment  is 
described  by  a  limited  number  of  layers  and  not  as  a  continuum. 

The  estimation  technique  needs  to  be  applied  to  a  larger 
number  of  data  sets.  By  examining  more  and  more  data  sets,  the  actual 
need  for  range  dependent  estimates  can  be  determined.  More  data  sets 
may  also  indicate  how  confidence  limits  can  be  defined. 
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